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We study the dynamics of the quantum phase distribution associated with the reduced density 
matrix of a system for a number of situations of practical importance, as the system evolves under 
the influence of its environment, interacting via a quantum nondemoliton type of coupling, such that 
there is decoherence without dissipation, as well as when it interacts via a dissipative interaction, 
resulting in decoherence as well as dissipation. The system is taken to be either a two-level atom 
(or equivalently, a spin-1/2 system) or a harmonic oscillator, and the environment is modeled as a 
bath of harmonic oscillators, starting out in a squeezed thermal state. The impact of the different 
environmental parameters on the dynamics of the quantum phase distribution for the system starting 
out in various initial states, is explicitly brought out. An interesting feature that emerges from our 
work is that the relationship between squeezing and temperature effects depends on the type of 
system-bath interaction. In the case of quantum nondemolition type of interaction, squeezing and 
temperature work in tandem, producing a diffusive effect on the phase distribution. In contrast, 
in case of a dissipative interaction, the influence of temperature can be counteracted by squeezing, 
which manifests as a resistence to randomization of phase. We make use of the phase distributions 
to bring out a notion of complementarity in atomic systems. We also study the dispersion of the 
phase using the phase distributions conditioned on particular initial states of the system. 

PACS numbers: 03.65.Yz, 42.50.Ct 



I. INTRODUCTION 

Open quantum systems are ubiquitous in the sense that any system can be thought of as being surrounded by its 
environment (reservoir or bath) which influences its dynamics. They provide a natural route for discussing damping 
and dephasing. One of the first testing grounds for open system ideas was in quantum optics [l[ . Its application to 
other areas gained momentum from the works of Caldeira and Leggett |2(, and Zurek [3|], among others. Depending 
upon the system-reservoir (S—R) interaction, open systems can be broadly classified into two categories, viz., quantum 
non-demolition (QND) or dissipative. A particular type of quantum nondemolition (QND) S — R interaction is given 
by a class of energy-preserving measurements in which dephasing occurs without damping the system. This may 
be achieved when the Hamiltonian H$ of the system commutes with the Hamiltonian H$r describing the system- 
reservoir interaction, i.e., Hsr is a constant of the motion generated by H$ [1, SUL A dissipative open system would 
be when H$ and H$ r do not commute resulting in dephasing along with damping 7] . A prototype of dissipative open 
quantum systems, having many applications, is the quantum Brownian motion of harmonic oscillators. This model 
was studied by Caldeira and Leggett [2| for the case where the system and its environment were initially separable. 
The above treatment of the quantum Brownian motion was generalized to the physically reasonable initial condition 
of a mixed state of the system and its environment by Hakim and Ambegaokar Q, Smith and Caldeira Q, Grabert, 
Schramm and Ingold (lCj . and for the case of a system in a Stern- Gerlach potential 11], and also for the quantum 
Brownian motion with nonlinear system-environment couplings among others. 

The interest in the relevance of open system ideas to quantum information has increased in recent times because 
of the impressive progress made on the experimental front in the manipulation of quantum states of matter towards 
quantum information processing and quantum communication. Myatt et al. Il3j and Turchette et al. [l4j have 
performed a series of experiments in which they induced decoherence and decay by coupling the atom (their system-.!?) 
to engineered reservoirs, in which the coupling to, and the state of, the environment are controllable. An experiment 
reported in Ref. fl5j demonstrated and completely characterized a QND scheme for making a nondeterministic 
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measurement of a single photon nondestructively using only linear optics and photo-detection of ancillary modes, to 
induce a strong nonlinearity at the single photon level. The dynamics of decoherence in continuous atom-optical QND 
measurements has been studied by Onofrio and Viola [l6| • In addition to its relevance in ultrasensitive measurements, 
a QND scheme provides a way to prepare quantum mechanical states which may otherwise be difficult to create, such 
as Fock states with a specific number of particles. It has been shown that the accuracy of atomic interferometry 
can be improved by using QND measurements of the atomic populations at the inputs to the interferometer [171 ]. 
QND systems have also been proposed for engineering quantum dynamical evolution of a system with the help of a 
quantum meter [l8j . In a recent study of QND open system Hamiltonians for two different models of the environment 
describable as baths of either oscillators or spins, an interesting connection was found between the energy-preserving 
QND Hamiltonians and the phase space area-preserving canonical transformations [l!| . 

A class of observables that may be measured repeatedly with arbitrary precision, with the influence of the mea- 
surement apparatus on the system being confined strictly to the conjugate observables, is called QND or back-action 
evasive observables [2(| [U [22|, 23 1 . Such a measurement scheme was originally introduced in the context of the 
detection of gravitational waves 24,[25|]. The energy preserving measurements, referred to above, form an important 



class of such a general QND measurement scheme. Since they describe dephasing without dissipation, a study of 
phase diffusion in such a situation is important from the context of a number of experimental situations. 

The quantum description of phases [26|,[I3] has a long history [28|, [29|, [3(J, HH [HJ . Pegg and Barnett [3lj . following 
Dirac ^28|, carried out a polar decomposition of the annihilation operator and defined a hermitian phase operator in 
a finite-dimensional Hilbert space. In their scheme, the expectation value of a function of the phase operator is first 
carried out in a finite-dimensional Hilbert space, and then the dimension is taken to the limit of infinity. However, it 
is not possible to interpret this expectation value as that of a function of a hermitian phase operator in an infinite- 
dimensional Hilbert space 1331. 34 . To circumvent this problem, the concept of phase distribution for the quantum 
phase has been introduced [33tl35|. In this scheme, one associates a phase distribution to a given state such that the 
average of a function of the phase operator in the state, computed with the phase distribution, reproduces the results 
of Pegg and Barnett. 

A study of the quantum phase diffusion in a number of QND systems was carried out in Ref. [36[ using the phase 
distribution approach. In this work we extend the above study to include the effect of dissipation on phase diffusion. 
Throughout this paper, the bath is assumed to be a collection of harmonic oscillators starting from a squeezed thermal 
initial state. An advantage of using a squeezed thermal bath is that the decay rate of quantum coherences can be 
suppressed leading to preservation of non-classical effects [33, [38|, [39| . It has also been shown to modify the evolution 
of the geometric phase of two-level atomic systems [4(|. The plan of the paper is as follows. In Section II, we recollect 
some results on the quantum phase distribution in QND systems from [36, 39]. We extend the previous expressions, for 
a single two-level atomic system, to the case of two two-level atoms and further plot the quantum phase distribution 
for ten two-level atoms. Following Agarwal and Singh (4l| we also introduce the number distribution and use it to 
discuss the complementary between the number and phase distributions. In Section III, we study the quantum phase 
distribution of a two-level atomic system interacting with its bath via a dissipative interaction. The evolution is 
governed by a Lindblad equation. The phase distribution is studied for the system initially (a) in an atomic coherent 
state and (b) in an atomic squeezed state. For the system in an atomic coherent state, complementarity between 
the number and phase distributions is discussed. In Section IV, the quantum phase distribution of the system of 
a harmonic oscillator, in a dissipative interaction with its bath, is obtained. In Section V, an application is made 
of the quantum phase distributions obtained for various initial system states and S — R interactions, to study the 
corresponding phase dispersion. In Section VI, we present our conclusions. 



II. QUANTUM PHASE DISTRIBUTION: QND 

Here we recapulate, from [361 ] , the results of Quantum Phase Distributions for a two-level atomic system as well as 
that of a harmonic oscillator which undergo interaction with their environments via a QND type of interaction. We 
consider the following Hamiltonian which models the interaction of a system with its environment, modeled as a bath 
of harmonic oscillators, via a QND type of coupling (39 |: 

H = Hs + Hr + Hsr 

= H s + Y,^ k blb k + H s J29k(bk + bl) + H 2 s Y,4^- (!) 

A; k k k 

Here Hs, Hr and Hsr stand for the Hamiltonians of the system, reservoir and system- reservoir interaction, re- 
spectively. The gfc's are dimensionless coupling constants. The last term on the right-hand side of Eq. (1) is a 
renormalization inducing 'counter term'. Since [Hs,Hsr] — 0, (1) is of QND type. Here Hs is a generic system 
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Hamiltonian which will be used subsequently to model different physical situations. Assuming separable initial con- 
ditions with the bath being initially in a squeezed thermal state and tracing over the bath degrees of freedom, the 
reduced density matrix of the system of interest S, in the system eigenbasis, is obtained as [3§| 



-%(E n -E m )t i(E 2 n -E 2 m ) v (t) 



x exp 



(E m - £„) 2 7 (t) P s nm (0), 



where 



and 



?(*) = -I^rf^sin(u, fe t), 



2 , , 2 



coth 



{e luJkt - 1) cosh(r fc ) + (e^ 4 - 1) sinh(r fc )e 



(2) 



(3) 



(4) 



For the reservoir R to be considered as a proper bath causing decoherence and (possibly) dissipation, we need to 
assume a 'quasi- continuous' bath spectrum with spectral density I{oj) such that for an arbitrary function f{uo) the 
continuum limit implies [6[ 



E §f(^) — J dul{w)f(u>). 



(5) 



We consider the case of an Ohmic bath with spectral density 

1(lu) — — ue ' c 

7T 



(6) 



where 70, having the dimension of 1/ (energy) 2 la], and oj c are two bath parameters characterizing the quantum noise. 
Using Eqs. © and © in Eq. ©, we obtain 



T)(t) = - — tan _1 (w c t). 
Using Eqs. ©, ^ in Eq. (JH) and using the T = limit, 7(f) is obtained as [3£ 



(7) 



j(t) = ^ cosh(2r) ln(l + Lj 2 c t 2 ) - — sinh(2r) In 

2-K 47T 



(l + 4^(1 -a) 2 ) 
{l + u 2 c {t-2a) 2 f 



To 

47T 



sinh(2r) ln(l + 4a 2 w 2 ), 



(8) 



where the resulting integrals are defined only for t > 2a [42j. Using Eqs. (O, ([6]) in Eq. ((4]) and using the high T 
limit, j(t) is obtained as [3!| 



7(i) 



lok B T 

lok B T 
2irTiu r 



cosh(2r 
sinh(2r) 



2w c itan 1 (oj c t) + In 



1 + UJ 2 c t 2 

Auj c (t — a) tan -1 (2u) c (t — a)) 



-4u> c (t — 2a) tan 1 (w c (t — 2a)) + 4aw c tan 1 (2aoj c ) 



In 



[l+w 2 (t-2a) 2 ; 
[l + 4cj 2 (i-a) 5 



hi 



1 



1 + 4a 2 uj 2 . 



(9) 



where, again, the resulting integrals are defined for t > 2a [421 ]. Here we have for simplicity taken the squeezed bath 
parameters as 



cosh(2r(oj)) = cosh(2r), sinh (2r(u>)) = sinh(2r) 
$(w) = au>, 



(10) 



where a is a constant depending upon the squeezed bath. Note that the results pertaining to a thermal bath can be 
obtained from the above equations by setting the squeezing parameters r and $ to zero. It is interesting to note that 
in the context of quantum information, the open system effect depicted in this Section can be modeled by a familiar 
quantum noisy channel, viz., the phase damping channel [40l l43l 
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A. Two-Level Atomic Systems 

Here we consider the case where our system S is a two-level atom. The system Hamiltonian Hs is 

huj , , 

Hs = —<Jz, (11) 

where a z is the usual Pauli matrix. The form of the system Hamiltonian Hs, Eq. (jlip . when substituted in Eq. ([1]) 
has been used in the context of quantum computation [i^, |4(| [4?J . In the context of a system of multiple two-level 
atoms, which is equivalent to an angular momentum system, we set Hs = TiloJ z . The Wigner-Dicke state (4§l. l49l [50j 
\j, m) , which are the simultaneous eigenstates of the angular momentum operators J 2 and J z , serve as the basis states 
for Hg 7 and we have 

Hs\j,m) = hujm\j,m) 

= E jim \j,m). (12) 

Here —j < m < j. Using this basis and the above equation in Eq. ((2]) we obtain the reduced density matrix of the 
system as 

Following Agarwal and Singh [4l[ we introduce the phase distribution V ((/>), <j) being related to the phase of the dipole 
moment of the system, as 

P(0) = ^±I/" d6sm{6)Q{6,4>), (14) 



4tt Jo 

where V{4>) > and is normalized to unity, i.e., J Q 27r d$V{<$) = 1. Here Q(9, </>) is defined as 

Q{e,<t>) = (0,0Ip s |0,0>, (15) 

where |0, 0) are the atomic coherent states [5l|, given by an expansion over the Wigner-Dicke states as 

\9,4>) = (j + m ) 2 (Bin(e/2)y +m (co S (0/2)y- m \j,m)e-^ +m ^. (16) 

m——j 

Using Eq. (|15p in Eq. (|14p . with insertions of partitions of unity in terms of the Wigner-Dicke states, we can write 
the phase distribution function as 



»M ■ 23 + 1 '" 



4tt 



/ d9 sin 6 ( e ^\ii n ) 

JO 



n,m——j 

x{j,n\p s (t)\j,m){j,m\9,<p). (17) 



Now we take up two physically interesting initial conditions for the system 5*. 



1. System initially in an atomic coherent state 

Here we consider the system S to be initially in an atomic coherent state which is the atomic analogue of the 
Glauber coherent state [49( . Thus the initial system density matrix is 

p s (0) = \a',[3')(a',f3'\. (18) 

Using the Eqs. (|13[) . (JTHJ) in Eq. (jTTJ) we obtain the phase distribution for a two-level atom, with j — i as [36] 



= tt - |~1 + t sin(a') cos(/3' +ut- fie-Wi® 

Z7T L 4 



(19) 



It can be easily checked that this V(<f) is normalized to unity. As can be seen from Eq. (|19p. only j(t) plays a role 
in the effect of the environment on the phase distribution. 
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2. System initially in an atomic squeezed state 



Now we consider our system S to be initially in an atomic squeezed state [51|, 152J, |53j, l54j expressed in terms of the 
Wigner-Dicke states as 

\(,p) = A p exp{<dJ z )exp(-i^J y ))\j,p), (20) 

where 

e 2e =tanh(2|C|) (21) 

and A p is usually obtained by normalization. Thus the initial density matrix of the system S is 

p s (0) = \C,P){C,P\- (22) 

Using the Eqs. (fT5|) . (I22p in Eq. (jTTJ) we obtain the phase distribution for a two-level atom, with j = i for p = ±^ 
as 36] 



Z7T 



1 ± . cos(^ - wtJe-^^W 

4cosh(9) v y 



(23) 



It can be seen that Eqs. (|23p. are normalized to unity. 

The above expressions may be extended to the case of multiple two- level atoms. For e.g., the quantum phase 
distribution for two two-level atoms, with j = 1, is: 

P(0 = JL jl ± — c 3 J h(2Q)) [cos(^ - art) cob(M V*)) cosh(e) 



sin 



(0 - wi) sin([^w] 2 ?7(t)) sinh(9)] exp(-[W] 2 7(^)) 



' cos(2[0 - tot]) exp(-4[H 2 7(*)) !> > (24) 



2(1 + cosh(29)) 



for p = ±1 and 



= £ {i - cos(2( ^ " «p(-4[H 9 7(*))} . (25) 

for p = 0. 

In comparison with Eq. (|23[) . which gives the quantum phase distribution for a single two- level atom, it can be 
seen that Eq. (|2"4"|) (phase distribution for two two-level atoms) involves both rj(t) and 7(f). This procedure may be 
carried to any number of two-level atoms using the Wigner-d function [55j : 

< p (t/2) = 2"V(i + n)!(i - n)!(j "TFT w ~ vj' ( 26 ^ 

where p (9) is the standard Wigner symbol for the rotation operator [55[ 

di p (9) = {j,n\e- iej «\j,p). (27) 

In Figure [1] we plot the quantum phase distribution for ten two-level atoms. It can be clearly seen from the figure 
that compared to the unitary case, interaction with the bath (characterized by finiteness of 70) causes phase diffusion. 
A comparison of the small- and large-dashed curves indicates that with increase in bath exposure duration t, the 
phase distribution diffuses as well as shifts to the right. It is also evident from the figure that increase in the bath 
squeezing r and temperature T also cause phase diffusion. The phase distributions are normalized to unity. 

In the case of QND type of interaction, the system is decohered without its energy being affected. This is reflected 
in the fact that with higher noise, the 'phase' gets completely randomized, resulting in a flattening of the distribution 
P(0), as depicted in Figure [T] whereas the 'number' distribution, given by 

p( m ) = (j,m\p s (t)\j,m), \m\<j, (28) 

. ^ ) (sin(a72)) 2 W+ m )(cos(a'/2)) 2 W- m ) for initial atomic coherent state Eq. fig]) 
J + m ' (29) 

\A p \ 2 e 2me \d : > mp (iT/2)\' 2 for initial atomic squeezed state Eq. 
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FIG. 1: Quantum phase distribution V{$) with respect to <j) (in radians) for ten atoms, starting in an atomic squeezed state 
(Eq. p2)n, with j = p = 5, a = (Eq. (TJOj)), 9 = -0.01832 (Eq. pi])! and 70 = 0.025, undergoing a QND system-bath 
interaction. Here ui = 1 and lj c = 100. The bold curve represents unitary evolution for t = 0.1, while the small-dashed and 
large-dashed curves are for the bath squeezing parameter r = 1.0, temperature (in units where % = ks = 1) T = 0.0 and 
evolution times t = 0.1 and 1, respectively. The dot-dashed curve represents the case r — 2.0, t = 0.1, T = 0.0, and the dotted 
curve the case r = 1.0, t = 0.1, T = 300.0. 



remains unaffected. The distributions p(m) and P{4>) may be thought of as complementary [4l| in the sense of 
conjugate Hermitian observables. For example, it may be verified that a 'number' state, i.e., Wigner-Dicke state, 
corresponds to a phase distribution of maximum uncertainty (in the entropic sense) [56(. This process may be 
understood as the selection of states in a preferred pointer basis [E, IE], which in this case are the Wigner-Dicke 
states, because of the nature of the system-reservoir interaction, whereby the environment 'monitors' the system in 
the preferred basis. As p(m) represents information in the preferred basis [E], the influence of the environment is not 
seen explicitly in Eq. (|29|) . 



B. Harmonic Oscillator System 



Here the system of interest S is taken to be a harmonic oscillator with the Hamiltonian 



// (a f a+-) (•'!()) 



The number states serve as an appropriate basis for the system Hamiltonian and the system energy eigenvalue (|30j) 
in this basis is 

£» = Mn + J). (31) 

Following Agarwal et al. HE] we define a phase distribution V{9) for a given density operator p associated with a 
state \9) as 

W = ^H >' O<0<2^, 



1 OO 

5- £ Pm.,ne l{n - m)b \ (32) 



2tt 

m.n— 



where the states \0) are the analogues of the Susskind-Glogower [23| phase operator and are defined in terms of the 
number states |n) as 

00 

= (33) 

n=0 

The sum in Eq. (I32p is assumed to converge and the phase distribution normalized to unity. Now we take up two 
physically interesting initial conditions for the system S. 
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1. System initially in a coherent state 
The initial density matrix of the system is 

P s (0) = \a)(a\, (34) 

where 

a = \a\e ie ° (35) 

is a coherent state [57J. Using Eqs. (|3"Tj) , (|34|) in Eq. ([2]) and then using it in Eq. (132)) . the phase distribution is 
obtained as 13611 



oo 

^ = 5 E 



1^1 pl (n- m )(e-e )-\ a \ 2 



III,!)- 



: yH!M! 



x ^-iu>{m — n)tgi(fiw) (m— n)(n+m+l)?7(t) e — (flu;) (n— m) 7(t) (36) 



2. System initially in a squeezed coherent state 
The initial density matrix of the system is 

p*(0) = |£,a)M, 
where the squeezed coherent state is defined as [13] 

|f,a)=£(Ol>(a)|0>. 



(37) 



(38) 



Here S denotes the standard squeezing operator and D denotes the standard displacement operator [57j| . Using Eqs. 
(l3"Tj) . (f3"7| in Eq. and then using it in Eq. (j3"2"|) . the phase distribution is obtained as [3(| 



i(n — m)8 



a i^(m—n) 



(tanh(n))" 



27r m tto~ 2^VH!(n!) C0Sh(ri) 



x exp [-|a| 2 (l - tanh(n) cos(20 o - -0))] 



xH„ 



ae 



i(9o-f) 



•y/sinh(2ri) 



tl„ 



ate 



y/ sinh(2ri) 



^,^y iuj(m n)tpi(huj) (m— n)(n+m+l)r;(i)g — (Ti.cj) (n— m) 7(4) 



(39) 



Here the system squeezing parameter £ = r\e 1 ^ and i? n [z] is a Hermite polynomial. The phase distributions depicted 
by Eqs. (|3"9")) have been plotted in Ref. f36{ . where they were seen to exhibit a phase diffusion pattern with the 

phase distributions being normalized to unity. 



III. QUANTUM PHASE DISTRIBUTION OF A TWO-LEVEL ATOMIC SYSTEM IN NON-QND 

INTERACTION WITH BATH 

Here we will obtain the quantum phase distribution of a two-level atomic system in an interaction with a squeezed 
thermal bath such that it undergoes both decoherence and dissipation. The reduced density matrix operator of the 
system S is given by @, [E3] 

+ 7o(iV + 1) (a_p s {t)a + - ^a+a-p s (t) - ^p s (t)a+^ 

+ 7Q N (a+p s (t)a- - ^<T_a + p s (t) - ^p s (t)<r^ 
- lQ Mcr + p s {t)<j + - lQ M*<j-p s {t)<j-. (40) 
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In the context of quantum information, the open system effect dep icted by Eq. (|40|) can be modeled by a familiar 
noisy channel called the generalized amplitude damping channel OIH, [44| for zero bath squeezing. For the case of 
finite bath squeezing and temperature, the corresponding noisy channel has been obtained by us recently [58[ and 
could appropriately be called the squeezed generalized amplitude damping channel. 

In Eq. ([40]) . 70 , having the dimension of (time) -1 , is the spontaneous emission rate given by 

^ = -isbr- (41) 

and <r_|_, cr_ are the standard raising and lowering operators, respectively given by 



CT + = |1)(0| = 7; {(?x+iOy) , 

= \0)(l\ = ~(* x -ia y ), (42) 



with a z being the standard Pauli operator related to the raising and lowering operators as a~\ = o z . In the above 
equations, [a, b] — ab — ba. In Eq. (|40|) 



and 



N = 7V th (cosh^(r) + s inh J (r)) + sinli (r), (43) 
M = -- sinh(2r)e i<E, (27V th + 1) = Re 1 *, (44) 

Nt h = • (45) 

e k B T _ I 

Here iV" t h is the Planck distribution giving the number of thermal photons at the frequency uj and r, $ are squeezing 
parameters. The analogous case of a thermal bath without squeezing can be obtained from the above expressions by 
setting these squeezing parameters to zero. Eq. (|4"P|) can be solved using the Bloch vector formalism (cf. 0], [40j| ) . 
However, the solutions obtained thus are not amenable to treatment of the quantum phase distribution by use of Eq. 
(Tr7|) . For this purpose we briefly detail the solution of Eq. (|4"0)l in an operator form. We closely follow the derivation 
given by Nakazato et al. [5!| and extend it to the case of a squeezed thermal bath. 
The Eq. (|4T))) can be written as 



dt 

where 
and 



^-p s (t) = Ap s {t) + p s (t)A^ + [ 1+ a-p s (t)<T + + 7 _a + p s (t)a. - lQ Mcr + p s (t)a + - 7 oMV_p s (t)a_] , (46) 



7+=7o(iV + l), 7- =70^, (47) 

A = -i 7 /3 -i( 7 + 2zw)a z , 

I = 7+ +7- =7o(2iV+l), 
7 = 7+ -7- = 7o- (48) 
The following transformation is now introduced in Eq. (|46[) : 



p s (t)=e At p I (t)e AH , (49) 



^(t) = 7 + a_p J (*)a + e-^ + 7-a+p J (*k-e 7< " IvMo+p 1 (1)^+^ - 7o M VV(^-e~ l2wt . (50) 
at 



yielding 



The solution of Eq. (|50|) is facilitated by the introduction of superoperators having the following action: 

V-P = £T-/9<7+, V+p = <T + p(T-, 

V a _p = a-pa-, V+p = <7+pa + . (51) 
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Using Eqs. (|5Tj) , Eq. (|50|) can be written as 



Integrating we get 

p T (t) = p^(0) + ^[ 7+ (e^*-e-^ t )P_ + ( 7+ e^ t +7-e-^ t - 7 ' 5 )P_P + ]/(0) 
+ 4j [7-(e 7+ * - e-T- 4 )^ + ( 7 _eT+ t + 7+ e^-* - >f)V + V-] /(O) 



(52) 



7o M 
7o M* 



sinh(crf) ^ t 
a 

sinh(ai) 



— !— (e lut {cosh{at) - — sinh(crf)} - 1)V1V! 
7nM a 



7o M 



' (e- lwt {cosh(at) + — sinh(ai)} - 



(53) 



where 



All the other terms are as given above. Using Eq. (|49| in Eq. (f53|) we finally obtain the solution of Eq. (|40|) as 



(54) 



= -p s (0)(l + e- f * + 2cosh(at)e-V) 

+ iCT 2 p s (0)cr 2 (l + e-^ t -2cosh(at)e- 3 ^ 1 ) 



sinh(at)e 2 



^P s (0) ( ^(l-e-^ t ) + —smh(at)e- x 5 
• ~ ' a 



7+ 



7- 



(1 - e"^) -^a_^(0)a+ + ^a + p s (0)a 



(55) 



This is the desired form of solution of the master equation (|4"01) . For the case of a thermal bath without squeezing, r 
and $ are zero and it can be seen that Eq. ([55)1 reduces to the solution obtained by Nakazato et al. [5^| for the case 
of a two-level atom interacting with a thermal bath. We will use Eq. ([55)) in the following subsections to investigate 
the quantum phase distribution. 



A. System initially in an atomic coherent state 

Taking the intial density matrix of the system S to be as in Eq. (TT5)) . using it in Eq. (f55j) . and then in Eq. (|17[) . 
with j = |, we obtain the quantum phase distribution as 

1 r 7r r 

V{4>) = — 1 + — sinfa )< acosh(ai) cos(4> — (3 ) + wsinhfat) sin(c/> — /3 ) 

- 7o i?sinh(ai) cos($ + + 0) j< 



(56) 



Here R, $ come from Eq. (|44[) . and 7^, a are as in Eqs. (|4"5| and (|54")) , respectively. The Eq. (f56"|) can be seen to be 
normalized to unity. When 7o is set equal to zero, i.e., for the case where the effects of the bath are neglected, Eq. 
(156]) becomes 



V{<f>, 7o - 0) = — {1 + j sin(a') cos(/3' + wt - </>)}. 



(57) 
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FIG. 2: Quantum phase distribution V((f)) (Eq. (J56J ) with respect to (j> (in radians), for a two-level dissipative system initially 
in an atomic coherent state (|18|) . Here uj — 1.0, <E> = tt/8, a' = /3' = 7r/4, and 70 = 0.25. The bold and small-dashed curves 
correspond to temperature (in units where % = fcg = 1) T — 0, bath squeezing parameter r = 0, but with bath exposure times 
t = 0.1 and 1.5, respectively. The large-dashed and dot-dashed curves correspond to T = 300 and t = 0.1, but r — and 2, 
respectively. Comparing the last two curves, we note that, counterintuitively, squeezing resists diffusion. 



In the analogous case of the QND system-bath interaction, the phase distribution was given by Eq. (fTO)) which with 
the bath coupling parameter 70 set to zero, is easily seen to reduce to Eq. (|57|) . This is a nice consistency check for 
these equations. 

Figure[2]illustrates the combined effects of temperature, evolution time and bath squeezing (r, $) on quantum phase 
distribution. Comparison of the small- and large-dashed curves brings out the diffusive influence of temperature, 
while a comparison of the bold and small-dashed curves shows that the phase distribution shifts with increase in 
bath exposure time. On the other hand, a comparison between the large- and dot-dashed curves illustrates an 
interesting feature of squeezing in dissipative systems governed by Lindblad-type equations (|40| . in that squeezing 
tends to counteract the influence of temperature, which in this case manifests as rcsistence to randomization of phase. 
A similar behavior is observed in the joint effect of temperature and squeezing on the geometric phase of a qubit 
(two-level system) interacting dissipatively with its environment (4oT | . The normalization of the phase distribution is 
preserved. 

We plot in Figure [3] the function 



p{m=\/2,t) = <l/2|p s (*)|l/2) 



l_7o 

7 /3 



1 + ^ 



sin V/2) + ^(l 



cosV/2). 



(58) 



Figure [3] depicts an expected behavior of a two-level system subjected to a dissipative channel. In particular, for 
T = and r = 0, it drives the system towards a pure state (with m = —1/2) and thus behaves as a quantum deleter 
[43I ]. Correspondingly, the phase distribution P(<fi) tends to level out for large bath exposure time t, as seen in Figure 
3J This brings out nicely the complementarity between p(m) and P(tp) |4l|. It is to be noted that, in contrast to the 
QND case, here the Wigner-Dicke states are not the preferred basis, and hence the environmental effects manifest 
themselves in the function p(m) as seen in Eq. ([58]) . 



B. System initially in an atomic squeezed state 



Taking the intial density matrix of the system S to be as in Eq. (f22|) . using it in Eq. (f55j) . and then in Eq. (|17j) . 
with j = ~, we obtain the quantum phase distribution for p = ±^ as 



2n 



1 ± 



7T f UJ 7o^? 1 

— ; — -< cosh(crf) cos(0) -I sinh(at) sin(^>) smh(cvt) cos(</> + $) > 

$h(Q) l a a > 



4 cos 



(59) 



Here is as defined in Eq. (f2"Tj) and all the other terms are as given above. The Eqs. (|B"9"|) are easily seen to be 
normalized to unity. Also by setting 70 to zero in them, they are seen to reduce to the cases of 70 set to zero in their 
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FIG. 3: The distribution p(m = 1/2, t) (Eq. (|58[H for a two-level dissipative system starting in an atomic coherent state 
(Eq. (|18Jl . as a function of time for different environmental conditions. The bold curve corresponds to temperature T = 100, 
70 = 0.0025, r = <f> = 0, ui = 1, a' = /3' = 7r/4, illustrative of a system becoming maximally mixed with time. The large-dashed 
curve corresponds to T = 0,7o = 0.025, r = 0, and depicts quantum deletion [4j]. The small-dashed curve represents the case 
T — 0, 70 = 0.025, r = 1. Here time and temperature are in units where h = ks = 1. 
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FIG. 4: Quantum phase distribution P(</>) (Eq. ([56}) with respect to (/> (in radians), for a two-level dissipative system starting in 
an atomic coherent state (Eq. (|18[) . at various times with temperature (in units where Ti = ks = 1) T = and bath squeezing 
parameters r = $ = 0, 70 = 0.025, u> — l,a' = f3' — n/4. The large-dashed, small-dashed and bold curves correspond, 
respectively, to evolution times t — 250, 50 and 10. The large-dashed curve depicts the randomization of phase distribution 
at long times. Comparison of this figure with Figure [3] clearly brings out complementarity between the 'number' and 'phase' 
variables. In particular, comparision between the large-dashed curves in both Figures shows how as the state becomes increasing 
pure, tending to m = —1/2, with time, the corresponding complementary distribution P((j>) level outs. 

QND counterparts, Eqs. (|2"3")l , respectively. This serves as a consistency check for these equations. On comparing the 
above equations, for the quantum phase distributions, with the corresponding ones for the case of QND system-bath 
interaction, these are easily seen to be more complicated. This is a reflection of the fact that the phase distributions 
developed in this section are for a process that involves both dephasing as well as dissipation, in contrast to the QND 
case, which involves only dephasing. 

We plot in Figures [5] the quantum phase distributions V(9) for a two- level system starting in an atomic squeezed 
state (|22p . An interesting feature in Figures [5] is brought out by a comparison of the bold and large-dashed curves. 
Squeezing is seen to have the effect of resisting the diffusive effect of temperature on the phase. This is similar to the 
behavior seen in Figure [JJ and suggests that this is a generic property of squeezing in a dissipative interaction. A 
comparison of the small-dashed and bold curves brings out the diffusive effect of temperature on the phase distribution 
while a comparison between the small-dashed and dot-dashed curves shows that the distribution shifts with time. 
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(a) 

FIG. 5: Quantum phase distribution V{$) (Eq. (|59[1 with respect to cf> (in radians) for a two-level system starting in an atomic 
squeezed state (g2|. Here w = 1.0, $ = tt/8, 9 = -0.01832 and 70 = 0.025. Figure (a) refers to p = | and (b) to p = -±. In 
both figures, the large-dashed and bold curves correspond to temperature (in units where h = ks = 1) T = 300 and evolution 
time i = 0.1. The bath squeezing parameter r is, respectively, 0.5 and 0.0. The small-dashed and dot-dashed curves correspond 
to T — and r = 0.0, with time t being 0.1 and 1.5, respectively. 



The phase distribution normalization is preserved. 



IV. QUANTUM PHASE DISTRIBUTION OF A HARMONIC-OSCILLATOR SYSTEM IN NON-QND 

INTERACTION WITH BATH 

Here we will obtain the quantum phase distribution of a harmonic-oscillator system, H s = huj{a^a + i), in a 
dissipative interaction with a squeezed thermal bath. The reduced density matrix operator of the system S, in the 
interaction picture, is given by [Tj, 1571 ] 



d_ 

db 



P s (t) = l0 (N + 1) (ap s (tW - \^a P s (t) - \ P s {t)a^ 
l0 N (a^p s (t)a - \acJp s {t) - ^p s (t)aa^ 



+ 



+ 70M* (ap s (t)a - \{afp s {t) ~ \p s {t){af 



(60) 



In the above equation, TV, M are bath parameters which will be given below and 70 is a parameter which depends 
upon the system-bath coupling strength. The Eq. (ffJO")) can be solved using a variety of methods (cf. 0], [57|). 
However, the solutions obtained thus are not amenable to treatment of the quantum phase distribution by use of Eq. 
(152")) . For this purpose we again briefly detail the solution of Eq. ([50)1 in an operator form. We closely follow the 
derivation given by Lu et al. (60j . The following transformations are introduced [oij: 



where 



Using Eqs. (fBTj) we get 



p' s {t) = S^Csp s {t)S{C), a =&(t)aS(0, 
5(C) = < ^ 



a = cosh(|C|)a - — - sinh(|C|)a f - 



Using Eqs. ([51]), {B3]) in Eq. (BO}, we get 

d 



dt 



P s (t) 



aK+ + (IK- + {a + (3)K + ^ p s (t), 



(61) 
(62) 

(63) 
(64) 
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where 

a = 70 iVcosh(2|C|) +7ocosh 2 (|C|) - ^ sinh(2|C|)(AfC* + M*(), 

(3 = 7o iVco S h(2|C|) +7osinh 2 (|C|) - sinh(2|C|)(MC* + M*Q. (65) 
The parameters involved in the above equation need to satisfy the following consistency condition: 

coth(|C|) + |^A/*tanh(|C|) =2N + 1. (66) 

It can be seen that 

M = isinh(2r)(2iV th + l)e 1 *, 

N = 7V th (cosh 2 (r) + sinh 2 (r)) + sinh 2 (r), 

^th = -jl , C = re < *, (67) 

e k B T - 1 

satisfy Eq. (|66p . In Eq. (|64|) . K + , and Kq are superoperators satisfying 

K+p s = ap' s a\ K_p s = a*<p s a, K p' s = -^ap' s + pVa + p' s ). (68) 
These superoperators can be seen to satisfy: 

{K_,K + }p s =2K oP s , [K ,K ± ]p' s = ±K ± p s , (69) 

which coincides with the commutation relations of the su(l, 1) Lie algebra. This brings out the intimate connec- 
tion between the solutions of the master equation (f60|) and the generators of the su(l, 1) Lie algebra. Using the 
disentangling theorems of the su(l, 1) Lie algebra, Eq. |64|) can be solved to yield: 

p' S (t) = ^ e V-(t)K- e \n(y (t))K 0e y + {t)K + ^ S{ q^ (70) 

where 

yo(t) = 
y+(t) = 





"iO* 
2 


-f3e~ 






7o 


a(e~ 


-70* 


-1) 


(J3e~ 


■7o« 


— a) ' 


(3(e- 


-70* 


-1) 


(pe- 


-7o t 





M*) - (71, 

Using Eqs. (f70|) . ([6"Tj) . the solution of Eq. ([60|) can be written as 

p s (i) = 5(C){e 2 5 i e 2 '- (t)K -e ln(! ' o(t))Ko e a + ( * )x +^ t (C)p s (0)5(C)}5 t (C)- (72) 

This is the form of solution of the master equation which we will use for investigation of the quantum phase distribution. 
We will use a special initial state of the system, the squeezed coherent state, 

p s (0) = K,v){v,C\, (73) 

where 

|C,ry) = S(C)£>(77)|0>. (74) 

Here |0) is the vacuum state and D(rj) is the standard displacement operator. Substituting Eq. ([731 in Eq. I|72p. the 
solution of the Eq. ([6H)l starting from the initial state ([73")) . following Lu et al. [6(j, is obtained as 



*(t) = 1 C -B(t)\n(t)\ 2 ^ ( 0ft) \ 1 

k 



X 



£ (f) (*) v / npi(r(*))*-'«(*)) fc - p ic^(*),o<p^(*),ci, (75) 
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where 



and 



\(,rj(t),l) = S(0\m,l) = S(OD(fj(tW), 



(76) 



(1 - e-w*), »j(t) = ij- 



7o 



1 + £(*))' 



(77) 



where /3 is given by Eq. (gS]). In Eq. (7BJ), D(fj(t)) = e v(tW-n*(t)a and D(f](t))\l) is known as the generalized coherent 
state (GCS) [62|,[63[ and thus the state \(, f}(t),l) would be the generalized squeezed coherent state (GSCS) [63]- The 
GCS's were introduced by Roy and Singh j62[, where they demonstrated that the harmonic oscillator possesses an 
infinite string of coherent states. We see from Eqs. ([75]) and ([73]) that under the action of the master equation (1601) . 
which is of a Lindblad kind, a harmonic oscillator starting in a squeezed coherent state ends in a mixture that can be 
expressed as a sum over GSCS. Thus the above case can be thought of as a concrete physical realization of GSCS. 

This is an example of ultracoherence pertaining to master equations governing the Lindblad type of evolution such as 
Eq. (|60|) . Ultracoherence refers to the structure induced into the Fock space F(J~t) , over a finite or infinite dimensional 
Hilbert space 7i, by the action of all canonical transformations, both homogeneous (e.g., squeezing operation) and 
inhomogeneous (Weyl operators) [H,[65|]. Starting from the squeezed coherent state (|?3"|) of the harmonic oscillator, 
obtained by applying the canonical transformation U = S(()D(rj) (|T4[) to the vacuum state, and applying a canonical 
transformation (|61[) to the master equation (|60p. results in a mixture of ultracoherent states, which in this case is the 
GSCS. 

Making use of the Fock-space representation of GCS [62j 



|n, a(t)) = e 



1=0 



E 77 L l n - n (\a(t)\ 2 )[a(t)] l - n \l}, 



(78) 



where L l ~ n (x) is the generalized Laguerre polynomial, and substituting Eq. ([75]) in Eq. (|32|) . reverting back to 
the Schrodinger picture, we obtain the quantum phase distribution of a dissipative harmonic oscillator starting in a 
squeezed coherent state (|73l) as 



V{6) 



1 _ P -/3(*)lrK*)l 

e -\m\ 1 = ^ e -Mm-n)t e Kn-m)6 y G* m (()G v n (Q 

(1 +/?(£)) ^ u ' mVS; ' VSy 

V II m,n u.v.k 



2tt 



Pit) 



x Lr l (\m\ 2 )L*;~ p {\m\ 2 )- 



i\ P \ 



S* (+\\v-p+k-l (fi(t)) u ~ l + k 



(f)*(t)) 



(79) 



In the above equation, G m>n (Q = (m\S(Q\n) and is explicitly given, with C = ne 1 ^, as [63 

(-l) p {(2p)l(2m)V 



G2m,2p — 

Similarly G 2m +i,2p+i(C) is given by 

G2m+l,2p+l 



(p)!(m)! \ cosh(ri) 
/tanh(n)\ (m+p) 



F? 



exp (i(m — p)4>) 

1 1 

-P, -m; 



(-l)P /(2p + l)!(2n» + l)! 



(p)!(m)! V cosh d (n) 



2' (sinh(n)) 2 



exp (i(m — p)<f>) 



(80) 



~ P ' m '2' (sinh(n)) 2 



(81) 



As has been pointed out in [63j, G m n is nonzero only for either m,n both even or both odd. For convenience it is 
sometimes assumed that <f> is zero and z = T\ is real. Here T\ = r, due to the initial condition (|73p and is the 
Gauss hypergeometric function [66j]. 
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FIG. 6: A comparison of the quantum phase distributions V{6), for a harmonic oscillator system starting in a squeezed coherent 
state, for QND system-bath interaction (Eq. (|39[) ) with that for dissipative system-bath interaction (Eq. (|79[) ). The former 
(latter) is represented by the dashed (solid) curve. In both cases, temperature (in units where Ti = ks = 1) T — 0, the squeezing 
parameters r = n = 1, bath exposure time t = 0.1, 70 = 0.025, u) = 1 . In the former case, ip = and oj c = 100, while in the 
latter, $ = 0. 



In Figure El we make a comparison of the quantum phase distributions V(0) for a harmonic oscillator system 
starting in a squeezed coherent state for QND system-bath interaction (Eq. ([55)) ) with that for dissipative 

system-bath interaction (Eq. (|79p). A comparison of the distributions brings out the differing effects of the two types 
of system-bath interactions on them. The phase distributions are normalized. 



V. APPLICATIONS: PHASE DISPERSION 



From the perspective of experiments, a relevant quantity is the quantum phase fluctuation, which may be quantified 
by the variance a 2 = {<f> 2 ) — {4>) 2 ■ For example, Ref. [r^ presents measurement of phase variance on atomic populations 
using interferometry improved by QND measurements at the inputs to the interferometer. However, this measure 
of phase fluctuation has the drawback that it depends on the orig in of the phase integration. A measure of phase 
fluctuation that avoids this problem is the dispersion D [53, H3, Ha [6t| , 



D = 1 



d<pe- itf "P(cp) 



(82) 



In this section, as an application of the phase distribution formalism employed above, we study the phase dispersion 
D from these distributions. We also evaluated the variance for these distributions (not presented in this work), and 
found that in certain cases, there is in fact a qualitative difference of the behavior of these two quantities. 

Figure [7] depicts the behavior of dispersion D of of a ten two-level atomic system, starting from an atomic 
squeezed state, interacting with a squeezed thermal bath via a QND interaction, with respect to environmental 
squeezing parameter r (Eq. (fTU)) ). The dispersion is found to increase with temperature and squeezing, tending to 
the maximal value of 1, corresponding to the uniform distribution P(<f>) = l/2n. This indicates that for a QND 
type system-bath interaction both temperature and squeezing have a similar effect of causing diffusion of the phase. 
Increasing bath exposure time t also leads to the effect of leveling out D to 1. At T = 0, this leveling out takes a 
much longer time on account of the logarithmic dependence of j(t) (Eq. ([8])) on t, indicating a power-law decay. 

Figure [5] is analogous to Figure [7J except that the dispersion of 4> is plotted with respect to system squeezing 
parameter £ (Eq. (|2Tj) ). As 6 appearing in the expression for P(</>) has a logarithmic dependence on Q (|2Tj) , dispersion 
is insensitive to change in £ over the plotted range. However, as expected, the dispersion increases with temperature 
because of the diffusive effect of temperature on the phase distribution. 

Figure[5]illustrates the behavior of dispersion D of </> of a harmonic oscillator starting from a squeezed coherent state 
interacting with a squeezed thermal bath via a QND interaction, with respect to environmental squeezing parameter 
r (Eq. (jlOp ). The dispersion is found to increase with temperature and squeezing, tending to the maximal value 1. 
Here the large-dashed curve, which represents the case of unitary evolution, shows no variation with respect to change 
in environmental squeezing parameter r (Eq. (jlOp ). as expected. 
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FIG. 7: Dispersion as function of environmental squeezing parameter r, for ten two-level atomic systems starting in an atomic 
squeezed state (Eq. (|22[) ), at various temperatures for a QND system-environment interaction. Here a = 0.0 [Eq. (|10[) ], 
70 = 0.0025, 6 = -0.01832, t = 1.0, j =p = 5, uj = 1.0 and lu c = 100.0. The bold, small-dashed, large-dashed and dot-dashed 
curves correspond to temperatures T (in units where % = Hb = 1) 0, 50, 100 and 1000, respectively. 
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FIG. 8: Dispersion as function of system squeezing parameter for ten two-level atomic systems starting in an atomic squeezed 
state (Eq. (|22[) ). at various temperatures for a QND system-environment interaction. Here a — 0.0, 70 = 0.0025, t — 1.0, 
j = p = 5, u> — 1.0, uj c — 100.0. The logarithmic dependence of 6 on ( (Eq. (|2ip ) implies a low sensitivity of the phase 
distribution to The bold, dashed, dotted curves correspond to the temperatures (in units where h = fee = 1) T = 0.0, 50.0 
and 100.0, respectively. The dot-dashed curve represents unitary evolution (70 = 0). 

Figure [10] depicts the behavior of dispersion D of <j> of a two-level system starting in an atomic coherent state 
interacting with a squeezed thermal bath via a QND interaction, with respect to environmental squeezing parameter 
r (Eq. (|10p V As before, dispersion is found to level out with increase in temperature and squeezing, tending to the 
value 1, which corresponds to a uniform distribution. We note that the pattern in this Figure is quite similar to that 
in Figure El whereas the use of variance for the data of Figure [TU] produces a qualitatively different pattern. 

Figure QT] shows the behavior of dispersion D of (j> of a two-level system starting in an atomic coherent state 
interacting with a squeezed thermal bath via a dissipative interaction, with respect to environmental squeezing 
parameter r (Eq. (fit))) ). While in the case of QND system-bath interaction (Figures [7l [91 and [TO]) , the dispersion 
is symmetric about r = 0, it is not so in this case of dissipative interaction. Further, unlike in the case of QND 
interaction, here increase in absolute value of squeezing (r) can cause a decrease in the dispersion. This illustrates 
the counteractive influence of the bath squeezing on the thermal diffusion of the phase distribution. This opposing 
behavior of temperature and squeezing seems to be generic to dissipative systems [40| . With increase in time t, phase 
tends to become randomized, increasing dispersion at any given squeezing towards the maximal value of 1, indicative 
of the washing away of the non-stationary effects due to the squeezed bath [64[. From Figure [U we see that increasing 
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FIG. 9: Dispersion as function of environmental squeezing parameter r for a harmonic oscillator starting in a squeezed coherent 
state (Eq. (J3TJ) ) at various temperatures for a QND system-environment interaction. Here u) = 1, uj c — 100, \a\ 2 — 5, 
70 = 0.0025, and t = 0.1. Here the parameter a — 0, and the system squeezing parameters are n = 0.5 and ip = tt/4. The 
bold, small-dashed and dot-dashed curves correspond to temperatures (in units where h = ks — 1) T = 0, 100 and 1000, 
respectively. The large-dashed curve corresponds to unitary evolution (70 = 0). 
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FIG. 10: Dispersion as function of environmental squeezing parameter r for a two-level system starting in an atomic coherent 
state (Eq. fllSp l. at various temperatures for a QND system-environment interaction. Here a = 0.0, 70 = 0.0025, t = 1.0, 
ijj — 1.0, uj c = 100, and oi = f3' = 7r/4. The bold, small-dashed, large-dashed and dot-dashed curves correspond to temperatures 
(in units where h = ks = 1) T — 0, 50, 100 and 1000, respectively. 

the bath exposure time (t) tends to shift and level out the distribution pattern. For finite temperatures, the latter 
effect predominates, and one observes a steady leveling out with time, with dispersion D tending to 1. Interestingly, 
the use of variance in place of dispersion for the data in Figure [TT] results in a qualitatively different behavior. 



VI. CONCLUSIONS 



In this paper quantum phase distributions of a number of physically interesting systems, interacting with their 
environment via a QND or a dissipative type of coupling, are analyzed. The system has been taken to be either a 
two- level atom (or equivalently, a spin- 1/2 system) or a harmonic oscillator with the environment being modeled as 
a bath of harmonic oscillators, initially in a squeezed thermal state, from which the common thermal bath results 
may be easily extracted by setting the squeezing parameters to zero. The phase distributions are explicitly evaluated 
taking into account the effect of the different environmental parameters on the dynamics of the system starting from 
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FIG. 11: Dispersion as function of environmental squeezing parameter r for a two-level system starting in an atomic coherent 
state (Eq. (|18p ). at various temperatures for a dissipative system-environment interaction. Here 70 = 0.0025, t — 1.0, u = 1.0, 
uj c = 100.0, $ = 7r/8 [Eq. (144 1 ], and a' = /?' = 7r/4. The large-dashed, small-dashed, dot-dashed and bold curves correspond 
to temperatures T (in units where ?i = fcs = 1) 0, 100, 300 and 1000, respectively. 

various initial states. 

In Section II, we recalled previous work on phase distributions for QND systems [36[ of two-level atomic systems 
(Section IIA) for different initial conditions of the system, starting (1) in an atomic coherent state, and (2) in an 
atomic squeezed state; and also of a harmonic oscillator (Section IIB) with the oscillator starting initially in (1) 
a coherent state, and (2) a squeezed coherent state. In Section IIA, some of the above results were extended by 
considering the phase distribution for multiple two-level atoms. In particular we studied, in Figure Q] the effect of the 
environmental parameters on the distribution for ten atoms starting in an atomic squeezed state and undergoing a 
QND system-bath interaction. The increase in bath squeezing r and temperature T causes phase diffusion while the 
increase in the bath exposure time t, causes the phase distribution to diffuse as well as shift. The phase distributions 
are normalized. We also introduced the number distribution p(m), expectation of the reduced density matrix p s (t) 
in the Wigner-Dicke states \j,m). By regarding the variables m and <f> as the 'number' and 'phase' of the atomic 
system, the relationship between the distributions p(m) and P(</>) may be considered as expressing complementarity 
in an atomic context. 

In Section III, the reduced density matrix of a two-level system interacting with a squeezed thermal bath via a 
dissipative system-bath interaction, resulting in a Lindblad form of evolution, was obtained, which reduces to the one 
found by Nakazato et al. [59[ for the case of a thermal bath without squeezing. This solution was used to study the 
phase distribution for the system, starting (1) in an atomic coherent state, and (2) in an atomic squeezed state. The 
phase distribution curves preserve the normalization of the distribution. The phase distribution exhibit diffusion as 
well as shift with time, as seen from Figures [5] and [5j An interesting feature that emerges from our work is that the 
relationship between squeezing and temperature effects depend on the type of system-bath interaction. In the case 
of QND type interaction, squeezing and temperature work in tandem, and produce a diffusive effect on the phase 
distribution. In contrast, in case of a dissipative interaction, with the reduced system dynamics governed by a Lindblad 
equation (|40|) . squeezing tends to counteract the influence of temperature, manifesting as a resistence to randomization 
of phase. This was noted for example in a comparison between the large- and dot-dashed curves of Figure^ and also 
in comparison between the bold and large-dashed curves in Figures[5l A similar behavior is observed in the joint effect 
of temperature and squeezing on the geometric phase of a qubit (two-level system) interacting dissipatively with its 
environment (4p| . Complementarity between the variables m and (f>, by a comparison of the distributions p(m) and 
P(4>), was brought out in an interesting manner for the case of a dissipative system-environment interaction and seen 
from a comparison of Figure [4] with Figure [3] In Figure [3l for the case where temperature T — and bath squeezing 
parameter r = 0, the system tends to the pure state \j — 1/2, m = —1/2), as seen by the large-dashed curve. This 
corresponds to the action of a quantum deleter (43j by means of an amplitude damping channel [44| . Correspondingly 
the complementary distribution P{4>) is seen to level out (the large-dashed curve in Figured]), indicating complete 
randomization. 

In Section IV, the quantum phase distribution for a harmonic oscillator in a dissipative interaction with a squeezed 
thermal bath, with the system starting out in a squeezed coherent state, was obtained. An interesting fact that 
emerged was that under the action of the master equation (|60[) . which is of a Lindblad kind, a harmonic oscillator 
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starting in a squeezed coherent state ends in a mixture that can be expressed as a sum over GSCS. A comparison of 
this distribution with that of the analogous case for a QND system-bath interaction (Figure brings out the differing 
effects of the two types of system-bath interactions on the phase distribution. 

In Section [Vj as an application we studied the dispersion of phase using the phase distributions conditioned on 
particular initial states of the system. In the case of QND system-bath interaction, the profile of dispersion D is 
symmetric about r = 0, as seen from Figures [7J l9l and ITD1 In contrast, the profile of dispersion is not symmetric in 
the case of a dissipative interaction (Figure [TTj) . indicative of the greater complexity of the latter type of interaction. 
Dispersion is a measure of phase fluctuations. Since the phase distributions used here are obtained taking the effect 
of the environment into consideration, the dispersions calculated using them would set a realistic estimate on phase 
measurements in a number of experimental scenarios. 

We hope that the treatment of phase distributions developed here would be of interest both from a technical point 
of view as well as in the context of experimental situations. 
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